Noise-induced volatility of collective dynamics 
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Noise-induced volatility refers to a phenomenon of increased level of fluctuations in the collective 
dynamics of bistable units in the presence of a rapidly varying external signal, and intermediate 
noise levels. The archetypical signature of this phenomenon is that -beyond the increase in the level 
of fluctuations- the response of the system becomes uncorrelated with the external driving force, 
making it different from stochastic resonance. Numerical simulations and an analytical theory of a 
stochastic dynamical version of the Ising model on regular and random networks demonstrate the 
ubiquity and robustness of this phenomenon, which is argued to be a possible cause of excess volatil- 
ity in financial markets, of enhanced effective temperatures in a variety of out-of-equilibrium systems 
and of strong selective responses of immune systems of complex biological organisms. Extensive nu- 
merical simulations are compared with a mean-field theory for different network topologies. 

PACS numbers: 05.40.-a,05.45.Xt, 89.65.-s 



I. INTRODUCTION 

Noise has effects a priori unexpected on the organiza- 
tion of complex systems made of interacting elements, 
as shown by stochastic resonance (SR) coherence 
resonance Q, noise-induced phase transitions Q, noise- 
induced transport [3| and its game theoretical version, 
the Parrondo 's Paradox 5] . SR occurs in a system when 
a small applied (sub-threshold) periodic signal is ampli- 
fied by the addition of noise and the maximum of ampli- 
fication is found for intermediate noise strengths. More 
generally, SR refers to the situation where noise and non- 
linearity combine to increase the strength in the system 
response^ Among others, SR was shown to appear in 
optical and magnetic systems @, Q, thought to be 
relevant from the Earth clima te M and the dynamics of 
ice ag es flOl , to neurobiology [U [l2| and visual percep- 
tions [13]. Generally SR is studied in bistable systems, 
where the amplification of a subthreshold periodic signal 
is achieved through the synchronization of noise-induced 
inter- well hopping of the dynamic variable and the driv- 
ing signal. The signal is maximally amplified when the 
level of noise is such that the Kramers time, which is 
the intrinsic lifetime associated with the noise-induced 
transition between the two stable states, equals half of 
the period of the external forcing. The bistability of the 
dynamic variable can be given either explicitly [14j -for 
low dimensional systems- or emerge from the interac- 
tion of the many constituents as for the magnetization in 
the Ising model in the ferromagnetic case flol [l6| . How- 
ever, systems of many interacting constituents may de- 
part from the paradigmatic setting of SR, as there is no 
inter- well hopping for the macroscopic observable, and 
thus the bistability is only preserved at the microscopic 
level. 

The Ising model, driven by a periodic signal, has been 
extensively studied in the realm of the kinetic Ising model 
and dynamical phase transitions [l5l - [l9| . When a spa- 
tially extended Ising system -for temperatures below the 



Curie temperature- is forced by a weak periodic influ- 
ence, the magnetization performs dynamics around its 
nonzero equilibrium fixed point, resulting in a nonzero 
time-averaged magnetization. For a given temperature, 
by increasing either the field strength or the period of the 
signal, the system becomes able to hop between the two 
symmetric equilibrium fixed points inducing a zero aver- 
age magnetization. The nature of this transition, from 
a nonzero to a zero average magnetization, can be man- 
ifold and depends on the control parameter. For weak, 
subthreshold forcing strength and finite-size systems - 
where inter-well hopping is enabled by the fluctuations- 
SR is at the origin of the transition. In large enough sys- 
tems -such that finite-size fluctuations can be neglected- 
and for intermediate field strengths, they experience a 
dynamical phase transition through a nucleation process 
|16| . For stronger forcing, the transition is forced by the 
exogenous field, with no contribution from endogenous 
factors. 

In this paper, we investigate the behavior of systems 
composed of many interacting constituents under the in- 
fluence of a time-varying external forcing. The Ising 
model framework is used as a generic example of such 
systems. In contrast to the classical SR studies, where 
the period of the periodic forcing is of the order of the 
Kramers time, we are interested in much faster signals. 
Additionally, aperiodic signals are included to the exter- 
nal forcing, a setup much closer to real world examples, 
which will yield some surprising differences to the cases 
involving periodic forcing. 

We find, for periodic and aperiodic signals alike, that 
for intermediate values of the noise intensity, the system 
dynamics shows a maximum in amplitude [H, |20| . Inter- 
estingly, the phenomenon of increased amplitude, which 
consists of an amplification of the signal for a periodic 
forcing, morphs into an increase of the system- wide fluc- 
tuations, uncorrelated with the signal for an aperiodic 
forcing. We call this phenomenon "noise-induced volatil- 
ity" (NIV). 
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There are many examples of systems composed of 
many interacting units that are subjected to a rapidly 
varying -periodic or aperiodic- common forcing. A first 
example refers to the empirical observations of strong 
amplifications of thermal noise into effective renormal- 
ized temperatures by quenched heterogeneities in mate- 
rials [2l| . in organized flows in liquids [22[ and in gran- 
ular media near jamming [23j j . We argue that NIV also 
provides a conceptual framework to model the immune 
systems of complex biological organisms, viewed as multi- 
stable complexes, which switch their mode of operation 
under the influence of noisy perturbations by pathogens 
and other stress factors |24| - t2a |. 

Another important application of the proposed mech- 
anism of volatility amplification can be found in finan- 
cial markets. The phenomenon of "excess volatility" [27j 
constitutes one of the major unsolved puzzles in finan- 
cial economics and refers to the ubiquitous observations 
that financial prices fluctuate with much larger ampli- 
tudes than they should if they obeyed the fundamental 
valuation formula, linking the share price of a company to 
it's expected future dividends and discount factors |28| . 
The model described below can be applied to represent a 
market of interacting investors, where the external forc- 
ing represents the news, i.e. the publicly available in- 
formation about the traded assets, that investors use to 
update their estimates of the asset's fair value. In addi- 
tion to the phenomenon of the increased volatility com- 
pared to the news amplitude, our framework allows us 
to address two other well known phenomena of financial 
markets. Namely the fact that the news are poor pre- 
dictors of future price changes [29| and the phenomenon 
of clustered volatility, quantified by the slowly decaying 
temporal dependence of volatility |30j | . 

We document the phenomenon of noise-induced 
volatility by numerical and theoretical calculations on a 
stochastic dynamical version of the Ising model on fully 
connected, regular as well as random networks, in the 
presence of rapidly varying periodic and aperiodic signal. 
NIV also constitutes a new indicator for an approaching 
phase transition [31j . 

This paper is organized as follows. In the following sec- 
tion, we introduce the model studied and the measures 
chosen to quantify the phenomenon studied. In Section 
IIII1 we revisit the case where the system is driven by 
a periodic forcing, focusing on the case of fast signals, 
by means of Monte-Carlo simulations and by means of 
an analytical approach. In Section IIV1 we present the 
main contribution of the paper: the study of the system 
driven by an aperiodic forcing. In Section [Vj we go be- 
yond the fully connected case, focusing on different net- 
work topologies and on a paradigmatic example of this 
phenomenon: the excess volatility in financial markets. 
Finally, Section fVTl presents a discussion and conclusions 
of the obtained results. 



II. MODEL DESCRIPTION AND DIAGNOSTIC 
VARIABLES 

Consider a system composed of N interacting units 
that can be in one of two states: s = ±1. The units 
are updated sequentially, randomly chosen at each unit 
micro-time 5 = 1/N, i.e. N updates are equivalent to 
one time unit at the macroscopic level. The update of 
the state s, of a given unit i from t to t + 5 is given by 

Si{t + S) = sign 

The value Si(t + 8) is determined by three competing 
contributions: (i) a common external dynamic forcing 
term f(t) (force, pathogens abundance, news); (ii) an 
annealed unit-specific term £i(t) that we will call noise 
(thermal fluctuations or threshold, intrinsic susceptibil- 
ity of a unit immune system compartment, investor id- 
iosyncratic opinion or private information) ; (iii) an inter- 
action term between units controlled by the amplitude 
K(t) (elastic coupling, feedback loops between immune 
system elements, social impact). 

The system's behavior will be investigated under the 
influence of two different types of external signals. To 
relate to the existing literature, we will use a smooth pe- 
riodic signal, f p (t) = Asm(ut), with period 2-k/uj and 
strength A; this implies that, when averaged over time, 
the standard deviation of the signal is er/ p = A/y/2. 
Along this paper, we denote by signal amplitude the stan- 
dard deviation of the signal, cry. As a periodic signal 
is a rather stylized and artificial setup, we will, in a 
later Section, also analyze the response of such a sys- 
tem to a stochastic process. The simplest choice of a 
stochastic process with tunable characteristic timescale 
is the Ornstein-Uhlenbeck (OU) process, which has ex- 
ponentially decaying memory and is defined by df ap = 
—dfapdt + AdWt, with mean, strength A, inverse time- 
scale 9 > and Wt is a Wiener process with normalized 
variance and zero mean. The asymptotic solution of the 
OU-process is 

f av (t)=A [ e-^-^dWr, (2) 

J — oo 

which gives a signal amplitude oy = A/y/29. 

The noise term of each unit in Eq. ([TJ fol- 

lows an independent stochastic process, whose values are 
drawn from the cumulative distribution function G(0, D) 1 
with zero mean ((£i(t)) = 0) and variance D 2 . Thus, 
(&(*)£,(*')) = D 2 5(t-t')5 lJ . If f(t) = and G(0,d) cor- 
responds to a Logistic distribution, the dynamical rule 
of Eq. (Jl} is equivalent to the kinetic Ising model with 
Glauber dynamics (cf. appendix) where D 2 is related to 
the temperature. 

In the interaction term in Eq. ([T]), the matrix of weights 
LUij defines the network connectivity between units, both 
in topology and relative strength. We assume that the 



3 



interactions between units are governed by connections 
that evolve much slower than the dynamics of the whole 
system. This amounts to considering a static network 
with fixed normalized weights ^2 ■ uiij = 1 . The effective 
coupling strength is given by K(t), which may depend on 
time to reflect global softening-hardening in rupture pro- 
cesses, evolving physiological states of immune systems 
and changes of social cohesiveness and/or social influence 
in financial markets. 

The macroscopic dynamics of the system is captured 
by the instantaneous "magnetization" 



m{t) 



which fluctuates around its mean value 



Q = (m{t)) t = - I m(t)dt 



(3) 



(4) 



where T is the duration of the simulation. We study the 
normalized standard deviation 
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yj((m(t) - Qf) t 

VUWh 



(5) 



of m(t), describing the "volatility" of the system dynam- 
ics scaled by the signal amplitude, 07. The correlation 
between the input signal and the magnetization, defined 
by 



P(r) = 



((r(t + r)-Q)/(t)) t j 



(6) 



provides an additional insight on the level of synchroniza- 
tion between the external influence and the overall sys- 
tem dynamics at a lag of r. The lag where the correlation 
is maximal will be called optimal lag, r* = max T p(r). 

In the case of the periodic signal, a common measure in 
stochastic resonance research is the spectral amplification 
factor (SAF)[13], 



R = 



S u [m(t)] = S u [m(t)] 



(7) 



which is the ratio of the power spectrum density of the 
magnetization S^fm^t)] over the power spectrum density 



of the driving signal, S u [f p (t)] 
driving frequency u>. 



aj = A 2 /4, both at the 



III. PERIODIC SIGNAL 

A. Simulations results 

First, we consider an homogeneous, complete, network 
(wy = l/(N — 1)) and a constant coupling strength 
Kit) = k = 1. The results reported below are not sig- 
nificantly different for random graphs with large average 



connectivity or when the connections allow for an unbi- 
ased statistical sampling within the population. As pre- 
viously said, we set G to be a Gaussian distribution with 
standard deviation D, and zero mean. Even though the 
system loses it's equivalence to the kinetic Ising model 
with Glauber dynamics, all the qualitative properties of 
the system remain unchanged. Without external forcing 
(A = 0), the system experiences, as for the equilibrium 
Ising model, a continuous phase transition at D c ~ 0.80fc, 
separating the ordered phase with two stable fixed points 
at ±Q(D), from the disordered phase, with a single sta- 
ble fixed point at Q(D) — 0. For the equilibrium case 
(A = 0), the dependence of Q{D) as function of D is 
shown by the continuous line in Fig. QJb). 

In Fig. [lja), we plot the spectral amplification factor 
as a function of noise strength for signals with different 
periods. The symbols are obtained by simulations of the 
model with 10 6 units. We observe that, even for rela- 
tively small periods, an increase of amplification exceed- 
ing one order of magnitude is achieved for a broad range 
of intermediate value of noise, the hallmark of stochastic 
resonance. Fig. [ljb) shows the average magnetization, 
Q(D), which is the usual order parameter in the kinetic 
Ising model studies, for the same signals as in Fig. [Ha). 
In section 1111 Bl we will develop a theory which shows 
that the global dynamics can be assimilated to a mo- 
tion within a potential that exhibits a transition from a 
monostable to a bistable regime. For fast signals, Q(D) 
has a value close to the equilibrium fixed point (con- 
tinuous line), suggesting that the dynamics of m{t) can 
be well described as fluctuations around this stable (or 
metastable) point, i.e. m(t) performs intra-well dynam- 
ics. For slower signals, larger fluctuations around the 
equilibrium fixed point are observed, as Q(D) vanishes 
already for D < D Cl indicating symmetric oscillations 
around m(t) = 0, i.e. m(t) performs inter-well dynam- 
ics. By Dt(A, uj, N), we denote the threshold value of 
the noise strength at which the potential barrier between 
the two equilibrium fixed points become small enough, 
such that the system performs inter-well hopping and 
thus Q(D) goes to zero. The noise strength threshold 
Dt approaches D c as either uj or N are increased or A is 
decreased. In the opposite limits, it will tend to zero. 

Independent on the driving frequency, the maximum 
in the amplification is always found at Dt- For fast sig- 
nals where Dt ~ D c , this maximum is observed at the 
equilibrium phase transition. This happens in the pres- 
ence of two competing phenomena near the equilibrium 
phase transition: on the one hand, a divergence in the 
susceptibility -making the system very sensitive to small 
changes in the external influences-; on the other, critical 
slowing down, which inhibits the reaction of the system. 
For slow signals, where Dt < D c , together with a more 
abrupt vanishing of Q, a pronounced jump in the spec- 
tral amplification factor R, defined in Eq. ([7]) , is observed 
at D = Dt, where the response of the system is greatly 
increased by the transition from intra- to inter-well dy- 
namics. For Dt < D < D c , the amplification decreases 
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FIG. 1. For periodic signals: different measures as function 
of the standard deviation of the noise in Eq. fTJ, D, for weak 
periodic signals with amplitude <jj = 0.01. The period of the 
signal is in the legend of (b). Symbols are obtained for simu- 
lations with a system size N = 10 6 and the lines are the result 
of the linear approximation presented in the main text. The 
vertical thin dotted line indicates Dt for cu = 27r/512, where 
the transition from intra- to inter-well dynamics occurs. The 
vertical thick dotted line indicates the critical value of noise, 
D c , where the phase transition takes place. For D £]Dt,D c [, 
m{t) performs inter-well dynamics, (a) Spectral amplification 
factor, R, as defined by Eq. (0. (b) Average magnetization, 
Q. For D < D c , Q = indicates that the system performs 
inter-well dynamics, (c) The normalized standard deviation, 
Op. (d) Instantaneous correlation between m and /, p(0), 
defined by Eq. 1(5)1. (e) Correlation between m and / at the 
optimal lag. (f) Optimal lag, r*, normalized by the period of 
the driving force. 
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FIG. 2. For periodic signals with amplitude <jf p = 0.02/\/2 
and different values for the standard deviation of the noise, 
D, specified in the legend. The spectral amplification factor, 
R, as a function of the period, 2n/u. Symbols are obtained 
for simulations of a system with N = 9 X 10 4 and the lines 
represent the linear approximation given by Eq. (I15|l . Three 
regimes of R can be identified as a function of the period: 
1) increasing amplification with increasing period, 2) plateau 
with intra- well dynamics, 3) for D < D c , stark increase of 
amplification due to inter-well dynamics. Similar results are 
found for larger system sizes, where the third regime appears 
for larger periods as the finite size fluctuations are reduced. 



with D, as the position of the minimum (±mo(D)) ap- 
proaches for D approaching D c from the left. 

Fig. |Tfd) shows the dependence of the correlation at 
zero lag on the noise strength. A minimum of instan- 
taneous correlation is observed at the same values of D, 
where the maximum in R occurs. This result confirms the 
existence of a double peak of the non-normalized instan- 
taneous covariance, as was found by Leung et al.0,Q. 

The effect of the signal frequency on the system behav- 
ior is shown in Fig. [21 where we plot the amplification R 
as a function of the period of the signal for different val- 
ues of D. For D > D Cl where the macroscopic system 
dynamics is described by a monostable potential, the de- 
pendence is composed of two regimes. The first regime is 
where the amplitude of the oscillations of m(t) increases 
with the period as m(t) is pushed for longer durations 
into one direction, allowing for greater deviations from 
the origin. For larger signal periods, R reaches a plateau, 
which constitutes the second regime, where the diffusive 
motion of m(t) is confined by the potential. The same 
behavior is observed for D <C D Cl where the potential 
barrier between the two minima cannot be overcome by 
the system, restricting the dynamics to intra-well mo- 
tions. Finally, for D < D c , the dependence shows a 
transition into a third regime. If the potential barrier 
is not too high compared with the noise intensity and 
the finite-size fluctuations, the system is able to perform 
inter-well dynamics for large enough periods of the ex- 
ternal driving. These inter- well dynamics are observed as 
a second rapid increase in R. The period at which this 
transition happens is the double of the Kramers time. 
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B. Analytical approach 

In order to understand these results, we now develop 
a mean-field theory, which becomes exact in the thermo- 
dynamic limit and for weak signal amplitudes. As our 
system is composed of many interconnected units, we 
can rewrite Eq. ([1]) by replacing the interaction term by 
the global instantaneous magnetization and by explicitly 
writing down Eq. ([T]) in the form of 



Si(t + 5) = 



+1 if $(t) > -km(t)-f{t) 
-1 if < -km(t)-f(t) 



Averaging over multiple noise realizations, the expected 
value for the state of the i-th unit, at time t + S is thus 
given by 

(si(t + 5)} 6 = 1 - G( - km(t) - f{t)) 
-G(-km(t))-f(t)) 
= l-2G(-km(t)-f(t)), (8) 

where G{9) is the the cumulative distribution function of 
the noise term i.e. G{9) — J_ dO' g(9'). Summing 

over all the units, we get that the updated instantaneous 
magnetization is exactly 



m(t + S) = m(t) + — [ Si (t + S) - Sl (t)\ . (9) 

By identifying 1/N = S as dt, averaging over the com- 
plete population and taking the continuous limit, the dy- 
namics reads 



(^r) = <-*(* + dt))-m(t). 



(10) 



With m(i) = (dm(t)/dt) and substituting Eq. ((HJ into 
Eq. (HOI), we get 

m(t) = -m{t) + 1 - 2G( - km(t) - /(*)), (11) 

which constitutes a general closed form evolution equa- 
tion for the magnetization of the system. 

From the point of view of the mean-field limit, the noise 
£i (t) can be either quenched or annealed, as the complete 
noise is condensed into the last term in Eq. (TTTj) . For 
A = (no external driving), the stationary solution of 
Eq. (|11[) gives the dependence of the equilibrium fixed 
point mo(D) as the solution of the implicit equation 



m (D) = 1 - 2G(-km (D)). 



(12) 



This solution exhibits a supercritical pitchfork bifurca- 
tion as a function of D, as expected for an Ising-likc 
system, which is displayed by the continuous line in 
Fig. [ljb). The critical parameter is found equal to 
D c = ky/2/ir, when £,i(t) is drawn from a Gaussian dis- 
tribution. 



The emphasis of this paper is on the system's reaction 
to fast, sub-threshold (A <C 1) signals, so that inter- 
well dynamics can be neglected. Thus, a perturbation 
expansion m(t) — toq + mi(t) up to first order yields 



d 
dt 



mi(t) = -r](D) mi(t) + 4>{D) f(t) + 0(mf ), (13) 



where <f>(D) = 2g(—kmo), rj(D) = 1 — 2kg(—kmo) and 
g = dG/d^. The dependence of <f> and 77 as a function of 
the noise strength is displayed by the dash-dotted lines 
in Fig.[ljb). Based on Eq. fp~3|) . 4>{D) can be interpreted 
as the attenuation of the signal by the noise in the in- 
dividual constituents of the system as <fi(D) < 1/k for 
any D. The parameter r/(D) can be understood as the 
strength of the restoring force that tends to bring m(t) 
back to its equilibrium value, mo, after being driven away 
by the influence of f(t). The larger 77, the closer the dy- 
namics of m(t) will be to tuq and the shorter will be 
the memory of m{t). In the particular case where f{t) is 
constant, m\{t) approaches the fixed point / (f)(D)/r)(D). 
Since 4>(D) remains finite when D passes through D c , it 
is the vanishing of ri{D) at D = D c and its smallness 
in the vicinity of D c that is at the origin of the ampli- 
fied volatility. Based on Eq. (|13[) . we can now compute 
the approximate value of the different measures for the 
external signals and compare them to the simulations of 
the actual system. 

For the periodic forcing, the dynamics of the magneti- 
zation yields 

A cj) 

m p(t) = —5 n [— w cos(wi) + rj sin(wt)] + mo. (14) 

r\ A + ur 

Together with Eq. ([7]) , this gives a spectral amplification 
factor equal to 



4 / A<f> 
p A 2 \ v 2 +uj 2 



22,2 
or + r/ z 



I] 2 + 10 2 



(15) 



Fig. Ufa) shows that, for fast signals (where Dt — D c ), 
the value of R obtained from this approximation matches 
well with the simulation results. Deviations from the 
approximation appear for slower signals when Dt does 
not coincide with D c and non-linear effects cannot be 
neglected anymore. 

As it can be seen in Fig. [lja), the spectral amplifica- 
tion factor can reach values above 100, showing that this 
system, even without considering inter-well dynamics, is 
able to show remarkable reactions to a weak forcing. Two 
distinct amplification mechanisms of sub-threshold peri- 
odic signals can be identified by comparing the simulation 
with approximation results. The first mechanism, being 
present for finite and infinite systems, is the increase of 
the output amplitude by the decrease of the value of rj: 
by reducing the restoring force of m(t), such that it can 
be further displaced from mo, the oscillation amplitudes 
are increased. The second mechanism is the amplifica- 
tion through inter- well jumps, which is only present in 
finite systems, as a sub-threshold driving force cannot 
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overcome the potential barrier without the existence of a 
source of fluctuations, like finite-size effects. Note that in 
the thermodynamical limit, where the approximation is 
valid for any frequency, at D c where rj{D c ) = 0, it follows 
from Eq. (TTSj) that the fluctuations of m(t) for nonzero 
frequencies will always be finite. 

In addition to the spectral amplification factor, R, we 
also measure the normalized variance of m(t), a, which 
measures the volatility of the dynamics, independent of 
the exact shape of the power spectrum. This measure is 
convenient as it can be used for comparison with the ape- 
riodic signal, for which R is not defined. From Eq. 
the normalized variance of r(t) is 



~ 2 



UJ 



IF {mp{t))t = 



m p (t) 2 dt 



(16) 



The equivalence between the eqs. (Il5j) and (|16|) is due to 
the use of a linear response approximation in the macro- 
scopic dynamics of the systems, neglecting the response 
at higher order harmonics of the driving signal. As a 
consequence, the approximation of the spectral amplifi- 
cation factor, R p , is better fitted by the simulations than 
a p . 

From Fig. [He), we see that the mean- field approxima- 
tion matches well the values of a obtained by simula- 
tions for intermediate signal periods. For large periods, 
the inter-well dynamics destroys the match, and for fast 
signals (small periods), the finite-size fluctuations over- 
shadow the fluctuations induced by the signals. 

The correlation between m p (t) and f p (t) is given by 



Pp(t) 



r\ cos(o;r) + ui sin(ojr) 
sjrf + UJ 2 



and, for the optimal lag, we obtain 



arctan(^) 



(17) 



(18) 



which follows directly from Eq. (|14l) . The correlation at 
zero lag is shown is Fig. [IJd). As for the results of R, 
the simulation results are well captured by the mean- 
field approximation for high frequencies and deviate due 
to inter- well dynamics for lower frequencies. As was ob- 
served in (U, a dip in correlation is observed for inter- 
mediate values of the noise amplitude. This dip occurs 
at Dt, concomitant with the maximum in the amplifi- 
cation measured by R. This apparent contradiction can 
be explained by the results shown in Fig. [TJf ) , which 
plots the optimal lag between m(t) and f(t) normalized 
by the period of the signal. At Dt, m(t) and f(t) are 
maximally lagged, reducing the correlation at lag zero. 
On the other hand, the correlation at the optimal lag 
has a value close to one, explaining the high amplifica- 
tion of the signal. This behavior can also be found in 
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FIG. 3. Size dependence of the maximal correlation between 
the magnetization and the signal with 07 = 0.01 and different 
time-scales for several noise intensities: D = 0.7 (full lines), 
D — 0.8 (dotted lines), D = 0.9 (dashed lines). Characteristic 
time-scales of the signal are specified in the legends, (top) 
Periodic signals, where p(r*) converges to 1 for any value of 
D in the thermodynamic limit. The rate of convergence in- 
creases with the signal period, (bottom) In case of aperiodic 
signals, the value of p(r*) does not converge to 1 in the ther- 
modynamic limit and depends on D and 9. At D = D c , the 
maximal correlation, p(j*), is significantly reduces compared 
to D c . 



the approximation, although there the maximum of am- 
plification and minimum of instantaneous correlation is 
found at D c as our approximation neglects inter-well dy- 
namics. From Eq. (jTTJ) and Eq. p8[) it follows that, in the 
thermodynamic limit, there exists an optimal lag, r*, for 
which perfect correlation is achieved for any frequency 
and any noise strength, i.e. P p (t*) = 1. 

Perfect correlation is however not achieved for any fre- 
quency in finite systems as shown by the dependence of 
p(r*) in Fig. Hfe), where the deterioration of the cor- 
relation with increasing driving frequencies is observed. 
The origin of this effect lies in the finite-size fluctuations, 
which vanish in the thermodynamic limit, as Fig. [3Jtop) 
shows. p{r*) converges to 1 for infinite systems, at a rate 
of convergence depending on D and the driving period. 



IV. APERIODIC SIGNAL 

We now turn our attention to the case where the com- 
mon forcing is aperiodic. In this section, we will con- 
sider an external force, which is described by the OU- 
process introduced in Eq. @ . Fig. 2] illustrates the typ- 
ical dynamic behaviors of m(t) for different values of 
noise strengths D for a single realization of the driving 
force fp(t). For D > D c , the magnetization fluctuates 
around m {D) = 0, with increasing amplitudes as D ap- 
proaches D c . The fluctuation amplitude deceases again 
when m(t) performs intra-well dynamics for D < D c , 
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m(t) 




FIG. 4. For aperiodic signals: Time evolution of the magne- 
tization, m(t), for different noise intensities D obtained with 
the same realization of the driving force f(t) and £i(t). N — 
10 4 , (Jap = 0.04, 9 = 1.0, k = 1.0, D = 2.0, 1.0, 0.8 (smaller to 
larger amplitude of m(t)'s fluctuations) and D — 0.7 (bottom 
curve fluctuating around m(t) = —0.6). 



where Q(D) £ 0. 

We will compute the same observables as for the peri- 
odic signal and investigate the differences. The formal so- 
lution of the linearized version of the dynamics, eq. (1131) . 
is given by 



mi(t) = Acf> 



-v(t-r) 



-9(t- 



'dW T >dr, (19) 



which describes the dynamics around niQ. The normal- 
ized variance of m\(t) described by eq. (fT"9f is now given 
by 



S-2 



V (o + v) ' 



(20) 



In Fig. [5ja), we compare the normalized variance ob- 
tained by means of numerical simulations with this the- 
oretical result. We find a very good agreement between 
the two for fast signals, with the same deficits due to 
inter-well dynamics as for the periodic case with slower 
signals. However, by comparing a ap for the aperiodic 
signal with the periodic case, a p , we observe a major dif- 
ference. Whereas for the periodic case, the volatility of 
the dynamics shows a finite maximum value at D c , the 
volatility diverges if the system is driven by an aperiodic 
signal as T)(D C ) = 0. This divergence of the normalized 
volatility, a ap , is not to be understood as an explosion 
of the dynamics, as m(t) cannot exceed [— 1, +1], It re- 
flects the immensely amplified reaction to a weak exter- 
nal forcing, consistent with the diverging susceptibility 
in equilibrium phase transitions. The fact that this di- 
vergence is absent for a periodic forcing stems from the 
discreteness of the power spectrum of the input signal. 

The correlation between the forcing f ap (t) and the 
magnetization, m(t + r), is given by 



\Jif + 8r] 



(r) + 0) e 



(21) 



and is shown in Fig. [5Jb) for zero lag, together with the 
simulation results, which are found in good agreement. 



P(0) 



Pit*) 




FIG. 5. For aperiodic signals: different measures as function 
of the standard deviation of the noise in Eq. (JTJ) , D, for weak 
aperiodic signals with amplitude <r/ p = 0.01. The inverse of 
the signal's time-scale, 9, is given in the legend. Symbols are 
obtained for simulations with a system size N = 10 s and the 
lines are the result of the linear approximation presented in 
the main text, (a) Normalized standard deviation, a, measur- 
ing the volatility amplification, (b) Instantaneous correlation 
between m and /, p(0), defined by Eq. ([6]). (c) Correlation 
between m and / at the optimal lag. In contrast to the peri- 
odic signal, a system driven by an aperiodic signal is not able 
to follow the signal, even at the optimal lag, which is well 
described by the mean-field approximation, (d) Optimal lag, 
t*, normalized by the time-scale of the signal. 



The optimal lag for which p a p{T~) is maximum occurs at 



V - J] 

yielding a maximal correlation of 



Pa P (r* p ) =2^ 



(22) 



(23) 



Here, we find the second major difference between the 
periodic and aperiodic driving. For the periodic sig- 
nal, it is always possible to find a lag at which the cor- 
relation between the forcing and the system's response 
is perfect. For the aperiodic signal (see Fig. [5jc) and 
Fig. [^bottom)), on the other hand, even in the ther- 
modynamic limit, the dynamics of the system can be al- 
most unrelated to the forcing. Perfect correlation is only 
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reachable for very slow signals, i.e. lim p ap (T*) = 1. 

Given that the forcing -an OU-process- has a continu- 
ous power spectrum, and that the response of the system 
is frequency-dependent, the spectrum of the macroscopic 
dynamics is distorted when compared to the one of the 
forcing, which has the effect of decreasing the correla- 
tion. Indeed, the system is only able to follow the part 
of the signal spectrum with frequencies lower than r/(D), 
which describes the rate at which the system can effec- 
tively react to external stimuli. As with decreasing 9, 
the contribution of lower frequencies in the signal's spec- 
trum is higher, the correlation for fixed D increases with 
decreasing 9. 

For D w D Cl the volatility amplifies many times that 
of the driving signal f(t). Concomitantly, p vanishes for 
every value of the lag r, indicating that the volatility of 
the system is generated by an internal collective behav- 
ior. It is important to note that, even though the system 
dynamics are endogenously generated, they are initiated 
by an exogenous driving of the system. This is further 
confirmed by the good agreement between the approxi- 
mation and simulations for fast signals. It is the shadow 
of the diverging susceptibility together with the vanish- 
ing rate of the reaction of the equilibrium model at D c , 
which is responsible for the observed NIV phenomenon. 



V. EXTENSIONS OF THE PHENOMENON 
STUDIED 

A. Different Networks 

To show that the NIV-phenomenon, characterized by 
the increase of volatility and decrease of correlation to 
the aperiodic forcing, is robust with respect to the struc- 
ture of the network, Fig. [B] shows the normalized volatil- 
ity a and the maximum in correlation p(r*) as a func- 
tion of D (the standard deviation of the noise term in 
Eq. (|TJ|) for different networks. We consider a two- 
dimensional regular grid with Moore neighborhood and 
random small-world connections with varying concentra- 
tion p w . Changing p w from to 1 interpolates between 
the regular 2D lattice and the completely random net- 
work. For each p w , the peak in volatility is still con- 
comitant with the vanishing of p at the threshold value 
Dt (p w ) ■ The noise intensity threshold Dt (p w ) is increas- 
ing in p w , as larger global interconnection enhances the 
cooperative organization, and larger noise is needed to 
destroy the ferromagnetic state. 



B. Excess volatility in financial markets 

By the definition of our model given by Eq. (fTJ). it is 
clear that it is also interpretable as a model of opinion 
dynamics, where Si(t) is the opinion of agent i at time t, 
in line with the established literature on discrete choice 
[33| l . The external forcing f(t) can be seen as the news 
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FIG. 6. For an aperiodic signal: (top) the normalized volatil- 
ity <7 and (bottom) the correlation at optimal lag, p{r*), as 
function of the standard deviation of the noise, D, for dif- 
ferent small-world random connection concentrations p w of a 
two-dimensional regular grid with Moore neighborhood. For 
p w — 0, the network is a 2D regular grid with Moore neigh- 
borhood. For p w = f .0, the network is a random graph with 
an average degree, d = 4. The other system parameters are 
N = fO 6 , 07 = O.Of , e = 1, k = 1. 



that are common to all agents, the noise £i(t) contains 
the agents' private information and the coupling term 
represents the social interaction between agents. The 
dynamics of the global opinion is then given by m(t). 

When applied to the social system of financial markets, 
the agents are investors and Si(t) corresponds to their 
opinion on whether the asset is over- or under-priced and 
hence to their willingness to buy (+1) or to sell (-1). The 
global demand is the given by m(t), which impacts on the 
price as 



log[p(t+l)] =logb(t)] 



m(t + 1) 
A ' 



(24) 



Here A represents the liquidity depth of the market, 
which is assumed constant and m(t)/X is the financial 
return, r(t), from period t to period t + 1. This equation 
expresses a linear market impact of the demands, which 
is a commo n hy pothesis in stylized models of financial 
markets [34], l35| . The results below do not change quali- 
tatively for more general non-linear impact functions [36j j . 

To apply our model to the financial markets, we use 
the coupling strength, k, instead of D as the control pa- 
rameter. Rather than assuming a fixed coupling strength 
for investors, we propose that the impact of colleagues' 
opinions on a given investor may be slowly varying with 
time. This effect reflects the fact that, in times of 
greater uncertainty, investors tend to be more influence- 
able by their surrounding (37j . There are many varying 
sources of uncertainty that impact financial markets, in- 
cluding the economic and geopolitical climate and past 
stock market performance. In the spirit of Ref. (38[, all 
these factors are embodied into the notion that K(t) 
undergoes a slow random walk with i.i.d. increments 
K(t + St) - K(t) - N(0,a k ), which is confined in the 
interval [k — Afc; k + Ak]. This later constraint ensures 
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FIG. 7. (color online) Upper panel: Sample dynamics r(t) 
(black bursty line) when the coupling strength K(t) of the in- 
teractions between agents undergoes a confined random walk 
(green) in [k — Ak; k + Ak] with Afc = 0.5 and step size 
crjb = Afc/yoOOO. Lower-right panel: quickly vanishing (resp. 
long memory of) the auto-correlation of r(t) (thin lines) (resp. 
\r(t)\ (thick line)) for two values of k and of Ak. Lower-left 
panel: NIV resonance in the presence of a time varying K(t), 
with Afc = 0.1 (circles), 0.2 (squares), 0.5 (diamonds). The 
other parameters are N = 10 4 , a ap — 0.04, D = 1. 



that social imitation remains bounded. We could have 
used an Ornstein-Uhlenbeck process or any other such 
confining dynamics, without changing the crucial results 
presented below. More complex models of sophisticated 
investors involve the strategic adaptations of the traders' 
propensity to imitate to the reliability of their colleagues 
in recent outcomes I39l44ll. 



By the mechanism of sweeping of the coupling strength 
K{f) close to the critical coupling strength k c (for fixed 
noise strength D) [iH E3|, we expect and find transient 
burst of volatility in response to the aperiodic driving 
force f(t) with constant amplitude and time scale. Fig. [7] 
shows a typical realization, where the return r(t) exhibits 
transient bursts, associated with excursion of K (t) in the 
neighborhood of k c . The lower- left panel of Fig [7] shows 
the robustness of the NIV phenomenon as a function of 
the average coupling k: even with a fluctuating K(t), a 
large volatility peak appears for intermediate values of k. 
The lower-right panel shows very short-range correlations 
of r(t) but very long-range correlations of the financial 
volatility \r(t)\ (another equivalent proxy for volatility), 
very similar to empirical observations of financial returns 
30] . Such long persistence of the volatility can be traced 
back to the slow diffusive nature of K (t) in line with the 
investors' slowly changing trust in the economy. From the 
previous section also follows that during times of crisis 
and strong social interaction (fc close to k c ), the dynamics 
is generated mostly exogenously, well in line with the 
documented inability of news events to explain large price 
movements 1291. 



VI. CONCLUSIONS 

In this paper, we have investigated the behavior of a 
system composed of coupled bistable units under the in- 
fluence of a -rapidly varying- common exogenous forc- 
ing and independent noise sources. Independently of the 
shape of the driving force, intermediate noise strengths 
trigger a strong level of fluctuations of the macroscopic 
dynamics around the critical value separating the ordered 
from the disordered phases. For a periodic forcing, this 
peak corresponds to a pronounced amplification of the 
signal, with a strong correlation between the macroscopic 
dynamics and the driving force at the optimal lag, the 
paradigmatic signatures of stochastic resonance. 

When the driving force is aperiodic a similar peak ap- 
pears, but here the amplitude of the fluctuations exceeds 
by far those observed for periodic signals. Coincidental 
with the increase of fluctuations, the correlation between 
the driving force and the system dynamics is completely 
destroyed. This shows that even though these fluctu- 
ations are induced by the common forcing, the macro- 
scopic dynamics has an endogenous origin. This phe- 
nomenon of noise-induced volatility contrasts with that 
of stochastic resonance, with the major difference being 
that it is not the signal, but the fluctuations that are 
amplified. 

Moreover, this phenomenon of noise-induced volatility 
also constitutes a new indicator for the approaching of a 
phase transition [3l"| , and it applies to a broader range of 
real-world systems due to the more common setup given 
of a coupled system driving by an aperiodic forcing. 

As an example of a system where this phenomenon 
can be observed, we have proposed the social system of 
stock markets, in which we have been able to not only ex- 
plain the excess of volatility observed in stock prices, but 
also the apparent absence of correlation between news 
and price changes and the persistence of volatility dur- 
ing times of crises. 



Appendix: Equivalence to the kinetic Ising model 
with Glauber dynamics 

A popular update mechanism in the kinetic Ising model 
literature was introduced by Glauber (44|. In it, the prob- 
ability for a spin to flip is given by 



1 



Pm 



aflAEi + I 



(A.l) 



where AEi is the energy gained by the system through 
the spin-flip and = 1/kT. With s l = ±1 and E { = 
~ s i{^2j KijSj + /), which is the energy of the state s,, 
Eq. (jA.ip can be rewritten as 



1 



(A.2) 
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where A = J^j Kij s j + /• With the transition rate given 
by Eq. (|A.2|) , we can compute the probability of being in 
state Si at time t + S by 



p(si] t + 5) =p(sn t)p Si ^ Si +p(-Si\ t)p^ Si ^ Si (A.i) 

1 



p(^ t)(i 

p(-Si\ t) — 



e s t 2(SA + l 
1 



) + 



e - Si 2PA + i 

= (p(su t)+p(-Si; t)J 
1 



1 



-S;2/3A I i 



e - Si 2/9A + I 
= p(Si), (A.4) 



which is independent of time and gives us the probability 
of finding spin i in state Sj. Eq. (|A.4[) can be rewritten 



a.s 



i(t + S) = 



with 
with 



+1 with 
— 1 with 



Prob = {e- 2f3A + l)- 1 
Prob = (e 2 ^ A + I)" 1 

Prob = 1 - F{-A) 
Prob = F(-A) 



(A.S) 



where F{x) is the cumulative density function of a Lo- 
gistic distribution with zero mean and variance 7r 2 /12/3 2 . 

The model studied in this paper, defined by Eq. ((T|), 
can be rewritten as 



Si(t + S) 



if 
if 



m > - 

m < - 

1 with Prob 
1 with Prob 



1 - G(-A) 
G(-A) 



(A.6) 



with G(x) being the CDF of the probability density func- 
tion of with zero mean and variance D 2 . By direct 
comparison of Eq. (|A.5I) and Eq. (IA.6I) , one can see that 
the model defined by Eq. fTJ) is equivalent to the kinetic 
fsing model with Glauber dynamics if the distribution of 
the noise is chosen to be a Logistic distribution. 
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